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Abstract 

We analyze the moments of the isosinglet generalized parton distributions H, E, H, E of 
the nucleon in one-loop order of heavy-baryon chiral perturbation theory. We discuss in 
detail the construction of the operators in the effective theory that are required to obtain all 
corrections to a given order in the chiral power counting. The results will serve to improve 
the extrapolation of lattice results to the chiral limit. 



1 Introduction 



In recent years one has learned that many aspects of hadron structure can be described in the 
unifying framework of generalized parton distributions (GPDs). This framework allows one to combine 
information which comes from very different sources in an efficient and model-independent manner. 
The field was pioneered in [2j and has evolved to considerable complexity, reviewed for instance 
in [1J E3 [7] . As GPDs can be analyzed using standard operator product expansion techniques 
PQIB], their moments can be and have been calculated in lattice QCD Lattice calculations of well- 
measured quantities can be used to check the accuracy of the method, which may then be employed 
to evaluate quantities that are much harder to determine experimentally. This complementarity is 
especially valuable in the context of GPDs, because experimental measurements as e.g. in jlUj may 
not be sufficient to determine these functions of three kinematic variables in a model-independent 
way. Moreover, several moments of GPDs admit a physically intuitive interpretation in terms of the 
spatial and spin structure of hadrons, see e.g. [21 ITT1 IT21 IT5| . 

A notorious problem of lattice QCD is the need for various extrapolations from the actual sim- 
ulations with finite lattice spacing, finite volume and unphysically heavy quarks to the continuum, 
infinite volume and physical quark masses. Simple phenomenological fits are often still sufficient in 
view of the general size of uncertainties, but with increasing numerical precision more reliable meth- 
ods have to be applied. Chiral perturbation theory (ChPT) provides such a method [Hj. Describing 
the exact low-energy limit of QCD it predicts the functional form for the dependence of observables 
on the finite volume and the pion mass and also the finite lattice spacing |16j . At a given or- 
der in the expansion parameter, ChPT defines a number of low-energy constants which determine 
each of these limits. Some of these constants are typically known from independent sources, and the 
remaining ones have to be determined from fits to the lattice data. The task of ChPT is thus to 
provide the corresponding functional expressions for a sufficient number of observables. In this paper 
we contribute to this endeavor by analyzing the moments of the isoscalar nucleon GPDs H, E, H 
and E in one-loop order. 

The analysis of pion GPDs in ChPT has been performed in several papers |17[ I18[ I19j . In the case 
of the nucleon GPDs, the chiral corrections have been calculated for the lowest moments |2()l 1171 12 lj 
in the framework of heavy-baryon ChPT, which performs an expansion in the inverse nucleon mass 
1/M. Due to the kinematic limit taken in this scheme, the sum and difference of the incoming and 
outgoing nucleon momenta p^ and p 1 ^ are of different order in 1/M. As a consequence, the nth 
moment of a nucleon GPD contains terms up to nth order in the 1/M expansion. Given the rapidly 
growing number of low-energy constants in higher orders of ChPT, it has been assumed that the 
chiral corrections can only be calculated for the terms of lowest order in 1/M, i.e. for the form factors 
accompanied by the smallest number of vectors (p' — p)^. This would be a serious setback for the 
program sketched above. The aim of the present paper is to show that the situation is much better. In 
particular, we find that the corrections of order 0(m- w ) and O(m^) to all form factors parameterizing 
the moments of chiral-even isoscalar nucleon GPDs come from one-loop diagrams in ChPT and the 
corresponding higher-order tree-level insertions. 

This paper is organized as follows. In Section[2]we recall the relation between moments of nucleon 
GPDs and matrix elements of twist-two operators and rewrite it in a form suitable for the 1/M 
expansion. In Section |21 we discuss the construction of twist-two operators in heavy-baryon ChPT 
and give a general power-counting scheme for their contribution to a given nucleon matrix element. 
In Sections 21 and we identify the operators that contribute to moments of GPDs at lowest order in 
the chiral expansion and give the results of the corresponding loop calculations. We summarize our 
findings in Section H3 
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2 Generalized parton distributions in the nucleon 



The nucleon GPDs can be introduced as matrix elements of nonlocal operators. In this paper we 
limit ourselves to the chiral-even isoscalar quark GPDs, which are defined by 
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where a sum over u and d quark fields on the l.h.s. is understood, so that H = H u + H d etc. Here 
a is a light-like auxiliary vector, M is the nucleon mass, and we use the standard notations for the 
kinematical variables 
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As usual, Wilson lines between the quark fields are to be inserted in Q if one is not working in the 
light-cone gauge a^A^ = 0. The x-moments of the nucleon GPDs are related to the matrix elements 
of the local twist-two operators 
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where = \{D^ — D^) and S denotes the symmetrization of all uncontracted Lorentz indices and 
the subtraction of traces, e.g. S t^, = \it iiv + t Ufl ) — jg^yt^x for a tensor of rank two. It is convenient 
to contract all open Lorentz indices with the auxiliary vector a, 
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and in analogy for O. The matrix elements of the operators © can be parameterized as |1J |S] 
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The moments of the above GPDs are polynomials in £ 2 , 
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The restriction to even k in (JSJ) and © is a consequence of time reversal invariance. 

To calculate the chiral corrections to the nucleons form factors we shall use the formalism of 
heavy-baryon chiral perturbation theory, which treats the nucleon as an infinitely heavy particle and 
performs a corresponding non-relativistic expansion |22) . The evaluation of nucleon form factors in 
heavy-baryon ChPT is simplified if one works in the Breit frame |23| . It is defined by the condition 
P = 0, so that the incoming and outgoing nucleons have opposite spatial momenta p' = —p = A/2 
and the same energy p' = p = M7, where 



7 = V / 1-A 2 /4M 2 . (7) 

In the heavy-baryon formalism the baryon has a additional quantum number, the velocity v, which 
in the Breit frame is v = (1,0,0,0). The incoming and outgoing nucleon momenta are thus given by 
p = Mjv - A/2 and p' = Mjv + A/2. 

The Dirac algebra simplifies considerably in the heavy-baryon formulation. All Dirac bilinears 
can be expressed in terms of the velocity Vn and the spin operator 

Using that (vA) = (vS) = 0, one finds in particular 

u(p')u(p) = 7 u v (p') u v (p) , 
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where the spinors 
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are normalized as u v (p, s') u v (p, s) = 2M5 s i s . With © one obtains the following representation for 
the matrix elements © in the Breit frame: 
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with 
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The definition of the E n and E n is conventional but might be confusing as E n is not the nth moment 
of E(x, £, t) etc. We nevertheless use this notation, in order to make it easier to compare our results 
with those in the literature. Notice that according to (p^the terms with E n>k in (|12|) are only nonzero 
if k is even, whereas those with M n ^-i, E n ^-i and M n ^-i are only nonzero if k is odd. We will 
evaluate these form factors in heavy-baryon ChPT. It is straightforward to transform back to the 
original form factors using 
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3 Twist-two matrix elements in heavy-baryon ChPT 



Heavy-baryon ChPT combines the techniques of chiral perturbation theory and of heavy-quark ef- 
fective field theory |22| (for a detailed review see Ref. [M )- The effective Lagrangian describes the 
pion-nucleon interactions in the limit when m n ,q <C M, where q is a generic momentum. In this 
situation the velocity v of the nucleon is preserved in the process. One introduces the nucleon field 
with velocity v as |22| 

N{x) = e- lMovx (N v {x)+n v (x)) , (15) 

where Mq is the nucleon mass in the chiral limit. The fields N v (x), n v (x) respectively contain the 
large and small components of the nucleon field and satisfy fN v = N v , fn v = —n v . Their Fourier 
transform depends on the residual nucleon momentum, i.e. the original nucleon momentum minus 
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Mqv. Integrating out the field n v (x), one obtains an effective Lagrangian for the pion-nucleon system 
which involves the nucleon field N v (x) and pion field tt(x), 

<C C fT = £"K + C-KN ) (16) 

where 

C n = + 4 4 ) + . . . , C„ N = + £® + . . . (17) 

are expanded in powers of q. The explicit expressions for the lowest-order terms read [21] 

4 2 ) = ^ Tr (dpU&rf + + rf x j) , 
C^l = N v {i(vV)+g A {Su)}N v , 

+ (c 2 - ^-) M 2 + c 3 + (c 4 + [S", 51 n^j iV„ (18) 

with U = u 2 = ex]){iTT a T a / F}, the covariant derivative V M = + T^, and 

r M = 5 (^^u + «^ ut ) = 7j?2 ^ ^ V^" + 0(vr 4 ) , 

u„ = i (u + 8„u - ufyut) =-jd^ a T a + 0{^). (19) 

The trace Tr and the Pauli matrices r a refer to isospin space. As is done in current lattice QCD 
calculations, we assume isospin symmetry to be exact here, neglecting the difference between u- and 
d-quark masses. The leading-order parameters appearing in (|18[) are the pion decay constant F 
(normalized to F ~ 92 MeV) and the nucleon axial-vector coupling both taken in the chiral limit. 
The field x implements the explicit breaking of chiral symmetry by the quark masses, and in the 
isospin limit can be replaced by x ~^ m2 1> where m is the bare pion and 11 the unit matrix in isospin 

(2) 

space. Estimates of the low-energy constants c, in the second-order Lagrangian PlL, which are of 

order 1/M, can be found in j2E]- We note that £ijV induces corrections to the nucleon propagator, 
which we treat as insertions on a nucleon line. They read — % ((vl) 2 — / 2 ) / {2Mq) and Aic\m 2 , where 
I is the residual nucleon momentum, and are to be multiplied with a nucleon propagator i/(vl + iO) 
from C}L on either side. The pion-nucleon vertices following from CzL can be found in Appendix A 

of [21. 

In the following subsection we discuss how to construct the operators in the effective theory 
that match the twist-two quark operators ©• Nucleon matrix elements in the Breit frame are then 
obtained as [201 

(p'\0\p) =M 2 Z N u v (p')G (r',r) Uv (p), (20) 

with the spinors u v and normalization J\f given in (|1U|) and 1)11(1 . Here Go(r',r) is the truncated 
Green function for external heavy-baryon fields N v , N v and the operator O in the effective theory. 
The residual momenta of the incoming and outgoing nucleon are given by 

r = p-M v = wv- A/2, r' = p' - M v = wv + A/2 (21) 
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with 

A 2 

w = M(j-l)+5M = -— - 4am 2 + 0(q 3 ) , (22) 

where 5M = M — Mq is the nucleon mass shift. Finally, Zjy is the heavy-baryon field renormalization 
constant, 

z " - 1 - - S iog £ - 8m2 48( " ) + o(,3) ■ (23) 

(3) 

where d^sC/-*) ^ s a low-energy constant in the Lagrangian C~L. As explained in |22] the corresponding 
operator is required for renormalization but does not appear in physical matrix elements. The value 
of (iJs (A*) can therefore be chosen freely (with different choices resulting in different values for other 
low-energy constants), and in it was chosen such that it compensates the log(m 2 /// 2 ) term in (|23|) 
at the physical value of m. Since we are interested in the pion mass dependence of matrix elements, 
we must explicitly keep the logarithmic term in Z^. For further discussion we refer to Section [3.21 



3.1 Construction of effective operators 

We now discuss how to construct the isoscalar local twist-two operators in the effective theory that 
match the quark-gluon operators 0(a) defined in and (|IJ). The relevant operators in the effective 
theory can be divided into two groups: operators 0^ which contain only pion fields (and couple to 
the nucleon via pion loops) and operators O n N which are bilinear in the nucleon field. The matching 
of operators thus takes the form 

0(a) = 0,(a) + O nN (a) , 0(a) O vN (a) , (24) 

where we have taken into account that there is no isoscalar pion operator of negative parity (i.e. no 
0,r(a)). The pion isoscalar operators O n (a) have been analyzed in several papers and 
we shall simply use their results. 

Let us now list the building blocks for constructing the operators O n N(a) and 7T j\i(a), which we 
collectively denote by Q(a), omitting the subscript irN for ease of writing. They should be bilinear 
in the nucleon field and should be tensors that have n indices contracted with the auxiliary vector a 
according to (J3J). To build tensors we have the following objects with Lorentz indices at our disposal: 
the velocity vector v^, the spin vector S^, the derivative <9^, and the antisymmetric tensor £pj,u\p- 

We 

recall that any Dirac matrix structure can be reduced to an expression containing the spin operator 
Sp, and that the metric tensor g^ can be omitted in the construction because the twist-two operators 
are traceless. Using the identities 

{SftjSu} = -{v^Vv -g^ v ) , [S ll ,S v ]=ie^ v x p v x S p (25) 

we can impose that Sp should appear at most linearly, or quadratically as the commutator [Sfj,,S v ]. 
Concerning the derivative <9„, we find it useful to have it acting either on single nonlinear pion fields 
u, v) in the combinations or given in fi9|). or as a total derivative on the product of all fields, 
or in the antisymmetric form d ^ = \[d ^ — d ^) on the product of all fields to its right or to its left. 
This will make it easy to keep track of factors A„ in the corresponding matrix elements, which play 
a particular role as we shall see. To give operators with the correct chiral transformation behavior, 
the derivative d must appear in the covariant combination V M = d ^ + r^. The fields and derivatives 
used in our construction are then any number of u^, V M and x± = u'X u =t UX u between the nucleon 
fields N v and N v , and any number of total derivatives acting on the operator as a whole. In the 
sense of ()19|) we henceforth refer to d^, and as "derivatives". They have chiral dimension 1, 
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whereas x± nas chiral dimension 2 and will not appear at the order of the chiral expansion we limit 
ourselves to. 

We can decompose the pion-nucleon operators Q n (a) as 

n 

Q n (a) = M n ~ k ~ x (av) n - k Q n , k (a) , (26) 

fc=0 

where Q n ,kifl) = a ^i ■ ■ ■ a n k Qn 1 "' lik does not contain any factors (av). The k external vectors a in 
Qn,k{o) can be contracted only with derivatives d^, V M , u p and the spin vector S p , or with the 
antisymmetric tensor. There can be at most one factor (aS) as discussed after (f2"5|) . so that Q n ^k{o) 
has to contain at least k — 1 derivatives. We can hence write 1 

Qn,k = MQn,k-l + Qn,k,0 + Qn,k,l + • • • , (27) 

where the operator Q n< k,i has chiral dimension k + i. Note that due to parity the number of factors 
Sfj,, Ufj, and e pu \ p must be even for O and odd for O. We remark that the contraction of a with the 
e-tensor involves at least two derivatives, given that we chose to replace its simultaneous contraction 
with v\ and S p by [S^^Sy] using (|25|1. As a consequence, the antisymmetric tensor does not appear 
in the operators with lowest chiral dimension for a given k. 



3.2 Tree-level insertions 



At tree level, the matrix elements of the effective operators between two nucleon states are easy to 
calculate. Since u p and Tn involve at least one or two pion fields according to (|19|) . derivatives in 



the effective operators are to be replaced as <9 M 



0, and V A 



-iwvu with w given 



in 1)22(1 . Notice that, while generically the derivative V M counts as 0(q) in the chiral expansion, 
the kinematics of the external nucleon momenta forces wv p to be of order 0(q 2 ). As a result, the 
leading-order contributions of the operator Q n ^ to the form factors in (|12j) come from the terms with 

With (|26j) one readily obtains 



maximum number of factors and no factor wv^ 



<p'\O n , k (a)\p) = (a&) k ~ l u v (p') 



(aA)E^ k + [(aS),(SA)}Mi% 1 



u v {p) 



{p'\O n>k (a)\p) L =? (oA)*" 1 ^) 



2M{aS) 



+ 



(aA)(5A) 
2M 



M. 



(o) 



n,k— 1 



Uv(p) , 



(28) 



where the superscript on each form factor indicates the term of order 0(q°) in its chiral expansion. 
At this order, the form factors E nk and M n ^-\ of the vector GPD are related to the matrix element 
of the operator O n i~o, since the nucleon matrix element of the operator O n ^-i is zero at tree level. 
As explained above, this operator contains a factor (aS), which due to parity must be accompanied 
by the axial field u p and hence does not contribute to tree-level matrix elements without external 
pions. For the axial vector GPDs one finds that the form factor E n ^-i (M ni fc_i) receives its leading 
contribution from the operator O n ^-\ (O n ^ : i), given the required number of factors A M in l)28[) . 

Beyond leading order, tree- level insertions contribute to the form factors starting at order 0(q 2 ). 
Contributions proportional to A 2 are due to operators with d 2 or to a factor w from operators with 
V, or to the kinematic factors^ m Q12|) and N in (|20jl . Contributions proportional to m 2 are due 
to operators with \+ or with V and from the wave function renormalization constant in (|2()jl. 

1 Instead of M one could also use Mo or F in 12611 and 1271 . since all are of the same order in chiral power counting. 
We find powers of M most convenient, because they also appear in the form factor decompositions (1I2> . 
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In the results of the following sections we explicitly include the terms proportional to g\ in the 
expansion (f2*5|) of Zn, whereas contributions from d^ 8 are lumped into the coefficients describing the 
m? corrections due to tree-level insertions. 

3.3 Loop contributions 

Let us now consider a loop diagram with the insertion of the operator Q n {a). One easily finds that 
the term M n ~ k ~ x (av) n - k Q„ ifc (a) in the sum (|26j) can contribute to the form factors in (|12|) which are 
accompanied by at least n — k powers of (av), i.e. to E njm , M W)m _i, E n ^ m -\ and M„ im _i with m < k. 
Chiral counting determines which terms can contribute to a given order. Namely, the contribution of 
the operator Q nt k,i in a loop diagram has chiral dimension 

D Kl =4L + (k + i)+J2 dim V n (j) + dim V nN (j) - 2I„ - I N , (29) 
i=i i=i 

where L is the number of loops and {k + i) is the chiral dimension of the operator insertion. V n (j) 
and K-at(j) respectively denote the jth vertex from C w and C^n in the graph. N n and N n N are 
the corresponding total numbers of vertices, and /„- and Jjv are the numbers of pion and nucleon 
propagators. 2 Using the relation L = I w + In — N w — N-^n (see e.g. j2U) and the fact that for our 
specific diagrams In = N n N, we can rewrite this expression as a sum of positive terms, which makes 
it easy to identify the different contributions at a given order: 

D kii = 2L-l + k + (i + l) + Y, (dimK(j) - 2) + J] (dimKiv(i) - l) ■ (30) 

i=i 3=1 

For each vertex we can insert either the lowest or any higher order, i.e. dimT/^j) = 2,4,... and 
dim V-knU) = 1) 2, .... Note that a loop diagram with chiral dimension D^^ generates contributions to 
a nucleon matrix element of order 0(q d ) with d > D^^. This is on one hand because of the explicit 
factors M and Zjy in (|20j) . and on the other hand because the sum + r'^ = 2wv^ is of order 0(q 2 ) 
and thus one order higher than the generic power associated with residual nucleon momenta. 

The form factors enter a matrix element multiplied by factors (a A) or (SA) as given in ()12|) . 
Taking these into account, one finds that the chiral correction from Q n ^k,i to E nyTn and M ni?n _i has at 
least order D^ i — m, while for the form factors E n ^ m _\ and M„ ]m _i it has at least order D^ i — m + 1 
and Dk t i — m — 1, respectively. This is a main result of our paper and allows one to determine which 
operators need to be considered to calculate the corrections to a form factor to a given order in the 
chiral expansion. Because Dk,% contains a term k and because of the constraint k — m > 0, the number 
of loops and the order of the chiral Lagrangian required to calculate the lowest-order corrections for 
a given form factor do not grow with m. Instead, a growing number of factors accompanying a 
form factor in the nucleon matrix element requires a growing number of derivatives in the operator 

Qn,k ■ 

As an application of our general result we find that the form factors E n ^ m and M njm _i can receive 
corrections starting at order 

• 0(q) from one-loop diagrams with insertion of the operator Q n ^ m -\ and leading-order (LO) 
pion- nucleon vertices, 

2 Note that a nucleon propagator correction from a higher-order Lagrangian counts as one (nucleon- nucleon) vertex 
with two nucleon propagators on either side, see the discussion after 1191 . 
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Table 1: Four- vectors and their products appearing in the numerators of the loop graphs of Fig. ^ 
NN vertices (arising from nucleon propagator corrections) are not explicitly shown in the graphs. 



derivatives in operator insertion 




Am 






and wv/j, 




Ufj, 




vertices 


irNN at LO 


SI 




ttNN at NLO 


(vl)(Sl) ± (vl){SA) 




NN at NLO 


(vl) 2 -1 2 ±IA- A 2 /4 



• 0{q 2 ) from the one-loop diagrams with insertion of the operators Q n ,m,o and Q n ,m+l,-l an d LO 
pion-nucleon vertices, and from one-loop diagrams with insertion of the operator Q n ^ m -i and 
one next-to-leading order (NLO) pion-nucleon vertex or nucleon propagator correction. 

In turn, the form factor E nm -\ receives corrections starting at order 0(q 2 ) from one-loop diagrams 
with leading-order vertices and insertion of the operator Q n ,m,~i- For M njJTl _i the discussion of 
corrections up to order 0(q 2 ) is more involved and will be given in Section [4.11 

To conclude the discussion of power counting, we consider the contribution to the form factors 
E n>m and M n>m -i of loop graphs with the insertion of the pion operators O n (a), see (|24|1. Repeating 
the above argument and taking into account that now In = N^n — 1, one finds that such diagrams 
have chiral dimension 

At = 2L - 1 + dim C + Yl ( dim ^l?) - 2 ) + Y, (dim^JvC?) " l) • (31) 

Given that the leading operator 0™(a) contributing to O n {a) has the chiral dimension n, one finds 
that it can contribute to the form factors E nm and M n m —i starting at order 0(q n m+1 ). Note that 
because of charge conjugation invariance the isoscalar pion operators 0™(a) have even n and that 
due to time reversal invariance the form factors E n ^ and M n & vanish for odd k. Together with our 
power-counting formula one thus finds that E nn gets contributions from C"(a) starting at order 0{q) 
and M nn _2 starting at order 0(q 2 ). All other corrections from operators O n (a) to form factors E n k 
and M n ^ start at 0(q 3 ). 

Let us now take a closer look at the one-loop graphs with pion-nucleon operator insertions, which 
are shown in Fig.^ With our construction of operators explained in Section fc.ll we can readily analyze 
the origin of factors A„, whose number determines to which form factor a graph will contribute. Using 
(vA) = (vS) = and the form (|18|) of the LO and NLO pion-nucleon Lagrangian, we find that the 
numerators of the loop integrals are composed as specified in Table ^ The denominators of the pion 
and nucleon propagators respectively are (I 2 — m 2 + i0) and (Iv + w + iO), so that the loop integration 
turns tensors Z w . . . l^ Lj into tensors constructed from and g^ u . A factor A M that can be contracted 
with or S^ (i.e. is not contracted to A 2 ) can hence only originate from total derivatives <9 M in 
the operator insertion and from an NLO pion-nucleon vertex or nucleon propagator correction. We 
will see that this reduces considerably the number of operators contributing to the leading chiral 
corrections of nucleon GPDs. 
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WV — y I + wv 



A 

2 



/ + tut; + 4? 



A 

2 



A 

2 



f / + WV - § 



A 

2 



a be 

Figure 1: One-loop graphs with the insertion of a pion-nucleon operator O n N(a) or O n N(a), denoted 
by a black blob. Not shown is the analog of graph c with residual momentum / + wv + A/2 of the 
intermediate nucleon line. 



4 Chiral corrections up to order 0(q 2 ) 
4.1 Axial- vector operators 

Using the formalism developed in the previous section, we now evaluate the form factors up to relative 
order 0(q 2 ). Let us start by giving the operators in O nj k,i that have the maximum number of total 
derivatives contracted with or . It will turn out that these are required to produce the factors 
of (a A) and (5A) in the form factor decomposition (|12|) . With the constraints of parity invariance, 
we find 

0„, fc _i(o) = (iad^NviaS)^ + . . . , 
O nA1 {a) = c n , fc {iad) k (id^)N v S»N V + . . . , (32) 

where the . . . stand for operators with fewer total derivatives. One has E^' k _^ = b n ^/2 and M^l_ 1 = 
2c nj fc for the tree-level contributions at order O(q ). From the time-reversal constraints on the form 
factors it follows that the low-energy constants b n ^ and c n ^ are zero for even k. 

As derived in Section fc.^l the leading chiral corrections to E n ^_\ come from one-loop graphs with 
LO pion-nucleon vertices and the operator O n ^-\. Since this operator does not contain pion fields, 
one needs to calculate only graph a in Fig. ^ One finds 

} + ^ k m) rn" + E^t + 0^), (33) 

where the terms going with m? and t originate from tree-level insertions as discussed at the end of 
Section 13.21 Here and in the following we use the subtraction scheme of |14| for the loop graphs, 
subtracting 1/e + log(4-7r) + vp(2) for each 1/e pole in 4 — 2e dimensions. The renormalization scale 
is denoted by /i, and the fj, dependence of the logarithm in (|33|) cancels against the fi dependence of 

— (2 m) 

El), , which we have not displayed for simplicity. Note that the bare parameters m, F, g& can be 
replaced with their counterparts at the physical point within the precision of our result. Since the 
nonanalytic corrections in (|33fl are independent of the moment indices n and k, they apply to the 
entire nucleon GPD H(x,^,t), 

^+m 2 H ( - 2 ' m \x,0+tH ( - 2 ^{x,0+O(q 3 ). (34) 



En,k(t)=E^ k {l 



3m z g z A 
(4vrF) 5 



log 



m 



+ 1 



fr(x,e,t) = frW(x,o(i-5^| 



(4vrF) 5 



log 



m 



+ 1 
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Let us now consider the chiral corrections for M n &_i. It follows from (j!2j) that the relevant 
diagrams have to produce a factor (aA) fc (SA). By power counting, the form factor M n ^_\ could 
receive corrections of order O(q ) from diagrams with LO vertices and the operator insertion O n k,-i- 
Similarly, corrections of order 0(q) could come from the diagrams with LO vertices and insertion of 
O n) k+i-i or O n ,A;,0) an d from diagrams with insertion of O n ^-i and one NLO pion-nucleon vertex or 
nucleon propagator correction. One finds no operator in O n> kfl that has k or more partial derivatives 
contracted with or 5^, and the same holds of course for U) k-i- According to our discussion in 
Section 13.31 the graphs just discussed can thus produce at most k vectors A^ (not counting those 
appearing in A 2 ) and hence do not contribute to M n ^-\. At order 0(q 2 ) there is a number of 
possibilities: 

1. graphs with LO vertices and insertion of O nt k+2-i, C>n,fc+i,o or O n jn t \. The insertion of O n ,jfe+l,o 
does not produce a sufficient number of factors A M , whereas insertion of O n ^+2-i gives a factor 
(aA) k+1 (aS), which contributes to the form factor E Uj k+i but not to M n ^-i- A correction to 
M n £_i is obtained from insertion of the operator O n ^ \ given in ()32|) . which already provides the 
tree-level term of this form factor. Only the loop graph in Fig. is nonzero for this insertion, 
and the result is analogous to the one for the contribution of C n fe _i to E n ^_i. 

2. graphs with one NLO vertex or propagator correction and insertion of O n ^+i-i or O ni ko- 
Insertion of O n k,o does again not provide enough factors of A^, whereas graphs with O n ^ + i j 
give zero due to time reversal invariance. This can be seen by direct calculation, or by noting 
that M n> jt_i is only nonzero for odd k, whereas the coefficient & n ,fc+i is only nonzero for even k, 
as remarked below (1321) . 



3. graphs with insertion of O n ^-\ and (%) two loops with LO vertices, or (ii) one loop with two 
NLO pion-nucleon vertices or nucleon propagator corrections, or (Hi) one loop with one NNLO 
pion-nucleon vertex or nucleon propagator correction, or (iv) one loop with a pion propagator 
correction from . The operator insertion provides k — 1 factors of A^, so that two more 
factors must be provided by the vertices or propagator corrections (without being contracted 
to A 2 ). This is not possible in case (i), because the LO pion-nucleon vertices only involve pion 
momenta and the pion momenta in a two-loop graph can be parameterized such that they are 
independent of A (as in the one- loop graphs of Fig. Likewise, a pion propagator correction 
in case (iv) does not depend on A and can therefore not contribute. In cases (ii) and (Hi) one 
obtains nonzero contributions from the graph in Fig. The NNLO vertices and propagator 
corrections follow from the Lagrangian C ^ given in [22j . We find that the only term providing 
the two required factors of A^ is the irNN vertex generated by 

9A 



4M 2 



N v \ ( VS) (u V) + ( Vu) (5 V) \ N v . (35) 



Note that this vertex does not introduce a new low-energy constant, similarly to the term 

(2) 

proportional to qa in E-Llr, which generates the ttNN coupling at NLO. These terms arise from 
the 1/Mq expansion of the leading-order relativistic pion-nucleon Lagrangian iV (iy — Mq + 
\9Ails)N, see e.g. (21]. 



Putting everything together, we obtain 



M n Jt) = M 



(4ttF) 

+ M^ k m) m 2 + M^h + 0(q 3 ), (36) 



2 

m 

lOg + 1 



p(0) rn 2 g\ m 2 



n 



wv — y I + wv wv + -j wv — y wv + y 

a b 
Figure 2: One-loop graphs with the insertion of the pion operator C"(a), denoted by a black blob. 

where the terms going with m 2 and t are due to tree-level insertions. With (|33j) . and (|14|) one can 
write for the isoscalar quark GPD E(x,£,t) 



E (x , (, t) _ m {x , o {i - ^ [io g ™ + 1] } _ o ^ io g ™ 

+ m 2 £( 2 ' m )(x,£) + t£ (2 ' t} (z,6 + 0(9 3 )- (37) 
4.2 Vector operators 

The analysis of the vector operators proceeds along similar lines. The operator H} k o reads 

O„,fe,o(a) = b n , k (iad) k N v N v + c n , k (iad)^ 1 (»0 M ) N v [(Sa), S»}N V + . . . , (38) 

where the . . . denote operators with fewer total derivatives. One finds E^u = b n k and M^t-i = ~ c n,k 
for the leading-order tree-level insertions, which implies b n ^ = c n> fc+i = for odd k. 

According to (fT2*)) the graphs that give chiral corrections to E n ^ or M n ^-i must produce k 
factors of contracted with or S^. With the constraints from parity invariance one finds that 
the operator O n ^-\ does not contain terms which have k — 1 or more total derivatives contracted 
with or S 11 . With the results of Section 13.31 this implies that E n f. and M nk -i do not receive 
corrections from pion-nucleon operator insertions at order 0(q), and that corresponding corrections 
at order 0(q 2 ) can only come from the diagram in Fig. with LO vertices. One finds that the 
one-loop contribution to the form factor E n }. is canceled by the terms proportional to g\ in the wave 
function renormalization constant (|23|) . For the form factor M n ^ one obtains a correction 

(39) 



M^{l- gi log — \ 



We note that for n = 1, k = this implies a chiral logarithm for the isoscalar magnetic form factor 

GM,s(t), 



cut) = /4 0) {i - ^ log ^} + G%?m* + G$gt + o ( , 3 ) 



(40) 



where is the isoscalar magnetic moment of the nucleon in the chiral limit and where we have 
added analytic terms due to tree-level insertions. The form (|4()[l is consistent with the result of the 
relativistic calculation in |29| . 
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One finally has to evaluate corrections due to the diagrams in Fig. [21 with insertion of the pion 
operator C"(o), where n is even. We use the representation of this operator given in [TU]. 3 



n—2 



0^{a) = F 2 Y,a n ,k{iad) k Tr (aL) {2iad) n - k - 2 {aL) + (aR) {2iad) n - k - 2 (aR) 



(41) 



fc=0 
even 



with L M = W dfJJ and = Ud^U'. As discussed in Section 1531 the corrections to M n ^ start 
at order 0(q 2 ) for k = n — 2 and at higher order otherwise. They are due to diagrams with LO 
vertices, so that only the graph in Fig. [2K contributes. This is because the leading-order ttttNN 
vertex corresponds to an isovector transition of the nucleon, as follows from ((18)) and the expansion of 
r in (|19|) . Combining the result with the correction in (|39|) and adding analytic terms from tree- level 
insertions, we obtain 

Mn, k (t) = <1 {l " g^r log ^} + <V- 2 MP*® + <f W + M&t + 0(, 3 ) , (42) 
where k is even and 

m 2 {rj) 



n-2 



sal y- r 



3=0 
even 



r) 2 



m (rj) log 



*9\ 



n „i 

Z~ ] 3{3 - 1) A^lj J i dv rf~ 2 m 2 ( V ) log 



J=2 
even 



m 2 (r]) 

~1T 



with 



m 2 (rj) = m 2 — - (1 — rj 2 ) 



(43) 



(44) 



Here A^j^ is the chiral limit of the form factors A™ k (t) describing the moments of the pion isoscalar 
GPD, 



/ dxx^H^ix^t) = Y,(^) k Al, k {t) 

J-i tr n 



(45) 



k=0 
even 



The relation to the low-energy constants a n> k reads \T3\ 



(46) 



which implies 



Q"n,n— k 



j=k 
even 



(47) 



j=0 
even 



The corrections to E n ^ start at order 0(q) for k = n and at order 0(q 3 ) or higher otherwise. At one- 
loop order we obtain 0{q) corrections to E n ^ n from graphs involving only LO vertices. Corrections of 

3 Note that the normalization of the twist-two operators used here differs from that in |19| by a factor of 2. The 
coefficients d„ t k have the same normalization here and in |19|. 
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order 0(q 2 ) involve either graphs with one NLO vertex or propagator correction, or graphs with LO 
vertices and the subleading part wv of the residual nucleon momenta, see the discussion after (f^Of) . 
Our final result including analytic terms from tree-level insertions is 

E n , k (t) = E^ k + 5 ri M^\t) + E^(t)\ + E ( " k m) m 2 + E^h + 0( q 3 ), 



where the order 0(q) correction reads 



q 2 n_2 fl 

E n ^{t) = -M{2m 2 -t)^_^n- 3 -2 j_ x 



drj 



even 



M(2m A - t) 



2 ^ ^g 2 



8(4vrF) 2 



3=2 
even 



7^(1 — if 
m{rj) 

1 - rf 



drj 

l rn{rj) 



and the order 0{q 2 ) term is 



3m 2 g 2 A 
(4vrF) s 



n-2 



l0g ~ / ,«n,n-j-2 



3=0 
even 



+ 



(47TF) 1 



n-2 



3=0 
even 



dr] if(l — rj 2 ) 



9l 
32 



2t log 



m 2 (r]) 



+ i 



(t - 2m 2 ) 2 
m 2 {rj) 

m 2 (r}) 

r — — — 

^ 9 



+ M 



c\m log 



m 2 (rj) 



+ 1 



-c 2 m 2 (r?) log 



m 2 (ri) 



c 3 m 2 (r]) lo, 



o 2 2 2 n 

6m g A . m \ - • . 

lQ g^ E 2 j m 

^ 3=2 
even 



2(4vrF) 2 



7T(0) 

n,n—j 



6 



(4^)' 



+ M 



m 2 (j]) 



3=2 
even 



2/ | log — ^ + 1 



(t - 2m : 
m 2 (r]) 



2\2 1 



qm log 



+ 1 



C2m 2 (j\) lo; 



m 2 {r\) 



c 3 m 2 (r]) lo. 



m 2 (r/) 



+ 



(48) 



(49) 



(50) 



The integrals over 77 in (|43|) , (|49|) and (|5U|) are elementary, but we have not found a simple closed form 
of the result for general n. In the next section we give explicit results for the values and derivatives 
at t = of the form factors. 

Our result for the form factor M n ^{t) disagrees with [HOI, where it was taken for granted that 
the only operators which contribute at order 0(q 2 ) are those which already appear at tree-level in 
the same form factor. As our analysis shows, this holds indeed in many cases but not in all. For all 
other form factors our results agree with [3U] where comparable. 4 For n = 2 our results for the vector 
operators also agree with those of Belitsky and Ji |21|. 5 

4 Note that |-':iU| gives the correction to E„, n at order O(q) but not at order 0(q 2 ). 

5 When comparing results, one must take into account that 12 II uses MS renormalization, where for each pole in 4 — 2e 
dimensions one subtracts 1/e + log(47r) + whereas we use the scheme of [14] and subtract 1/e + log(47r) + tp(2). 
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5 Results for moments of GPDs 



We now transform the results of the previous section to the basis of the form factors A n ^, B n ^, C n 
and A Ut k, B n ^ corresponding to moments of GPDs in the conventional parameterization. We give 
the values and derivatives of these form factors at t = 0, which allows us to obtain closed expressions. 
Furthermore, these quantities are of most immediate interest in studies of GPDs on the lattice. 

With our results (|42[). (|ig|) and the conversion formulae f|13j). (jTl|) one obtains for the 

form factors at t = 



A n>h (0) 

B nt k(0) 

A n , k (0) 
C n (0) 



A 



B 



(0) 
ri.k 

(0) 
n,k 

(0) 



3m 2 g\ 
(4ttF) 2 
3m 2 ^ 
(4vrF) 2 

m 2 + 0(m 3 



1 



log 



A«J , /i(2,m.) 



+ 1 



m 

log — + 1 



B 



(0) 
n,fc 



«(0) 

71, fc 



(o)\ 3m 2 g 2 m 



+ 4 2 , m) m 2 +0(m 3 



7(0) ^ 2 ffj 

n,fe (A^T?\2 



m 



(lrrf;2 log-^ + ^V + O^ 3 ) 



(4vrF) 2 



log ~~2 + k,n-2 M^(0) + B^frr? + 0(m 



+ ^ 1,7r) (0) + 4 2 '" } (0) + C^m 2 + 0(m 3 ) 



(51) 



m 

J n,k' 



j(0) 



M 



(0) 



lp(0) 4 (0) 



(0) 



with coefficients related to those in Section 0] by A® k - ^ n fc , ^ n k - ^ n k - ^ n>k , - n>fe - - n , fc , 

E?^ ], = M^ ], — E^,, C^ ' = E^ and by analogous relations for the coefficients with super- 
script (2,m). Setting m, gA, F to their physical values and choosing /i = M, one finds that the 
corrections from loop graphs with nucleon operator insertions in (|51|) are moderately large, with 
3m 2 g A (4vrF)- 2 [logfm 2 /^ 2 ) + 1] « -0.20 and m 2 ^ (4vrF)~ 2 log(77i 2 //j 2 ) -0.09. In the case of 
E? n ^ this loop correction can be substantial if |-B n fc| <C |^4n,fc|> which is empirically found for the 
electromagnetic form factors (i.e. the case n = 1) and also in lattice evaluations (Hj for the moments 
with n = 2. The contributions to E? nin _2(0) and C n (0) from loop graphs with pion operator insertions 
are 



E^(0) 
E^(0) 



6m 2 g\ 

(4vrF) 2 *~ 6 /i 2 



J=2 
even 



3nmMg A 
2(4vrF) 2 



i=2 

even 



i + i 



7T(0) 

n,n—j ' 



3m 2 g% 
2(4vrF) 2 



+ 



12m 2 

(47TF) 2 



4 



M 



i=2 

even 



771 



Cl fog — o- + 1 



3 771 
-C 2 log — 

4 /i' 



C3 lo; 



711 



/i 



2 + 2 



3=2 
even 



j + 1 



tt(O) 



n,n—j 



(52) 



Setting M, m, ^a, E 1 to their physical values, choosing /j = M, and taking the estimates c\ 



-0.9 GeV" 



C2 
n(2,7r), 



3.3 GeV" 



4.7 GeV -1 from [251 we find M, 



r(2,7r) 



(0) 



-0.27 A 



E^' 71 ' (0) + E^'^O) ~ (0.12 + 0.17) A^o^- At the physical point, the order 0{m) correction is hence 



7T(0) 

2,0 



and 
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not very large. The full size of the order 0(m 2 ) corrections depends of course on the analytic terms 
in (|51|) . whose values are not known. 

The derivatives of the form factors at t = obtain nonanalytic contributions only from the pion 
operator insertions. Writing dtA(0) = [Jj-<4(i)] t=0 etc. we have 



d t A nik (0) = E^f + 0(m) 



n,k 
F(2,i) 



d t B n , k (0) = M%> - E%$>/2 - ^/(32M 2 ) + 0(m) , 



(0) 



cM n ,fc(0) = 



.(2,*) 
n,k 



ill 



(0) 



<l)/(4M 2 )+0(m), 



n,k 



d t B n , k (0) = fc,„- 2 d t M^(0) + M™ - Effi + (M® - <i)/(4M 2 ) + 0( 
d t C n (0) = d t E^(0) + dtE^Ho) + Egg + 4°)/(4M 2 ) + O(m) 



(53) 



with 

d t M^\0) 
2 

+ 



3gj 

(4^F) 2 

M 7T5 2 



m 

log — + 1 



.4 



tt(O) 



J (5j + 14) 



i=2 

even 



m 8(4vrF) 2 ^_ 

even 



tt(O) 



4(4vrF) 2 



m 

log — o- + 3 



(i + i)(i + 3)" n ^' 



even 



(4vrF) 2 



^ + M 



3 / m 2 

n - T c 2 ( log — + 1 



j _|_ 1 



r :; | log^- 



E 2 ~ 

i=2 

even 



j (i + 4) 



tt(O) 



(j + l)(j + 3)"™'^ 
(54) 



Note that in the chiral limit the derivative <9i_B nin _2(0) ~ dtMn' n \o) diverges as log(m 2 /^ 2 ) and 



dtC n (0) ~ dtEn' n \o) as 1/m. With the parameters specified above, one finds dtM^' 71 ' (0) « 
1.7 GeV" 2 A^ ( 0) and a t ^ 1,7r) (0) + dtE^ v \o) « -(2.5 + 1.2) GeV" 2 ^ ( 0) . Numerically, the term 
dtEn' n \o) is thus important but not extremely large at the physical point. 



(2,7T), 



6 Summary 

Using heavy-baryon chiral perturbation theory, we have calculated the chiral corrections up to order 
0(q 2 ) for the form factors which parameterize moments of nucleon GPDs. We have restricted ourselves 
to vector and axial-vector quark distributions in the isosinglet combination. Our results generalize 
trivially to the corresponding gluon GPDs, which have the same quantum numbers and therefore the 
same corresponding operators in the effective theory (except for the values of the matching constants). 
Our method is also applicable to operators of different tensor or flavor structure. 

The moments of GPDs contain terms of different order in 1/M, ranging from M n_1 to M~ l . We 
have shown that, due to the way in which factors Vn and A^ arise in the calculation, the number of 
loops and the order in the expansion of the effective Lagrangian required to calculate a form factor to 
a given order 0(q d ) does not grow with the number of factors A^ that accompany the form factor in 



16 



the nucleon matrix element. A general power-counting formula is given after (j30fl . In the case of the 
form factors M n ^{t), calculation of the order 0(q 2 ) correction requires the pion- nucleon Lagrangian 
up to third order. 

We have found that the form factors E n ^(t) and M n ^(t) receive corrections of order 0(q 2 ) which, 
apart from analytic terms, are independent of the moment indices and independent of t. The same 
holds for the one-loop corrections to M n ^{t) with k < n — 2, whereas the corresponding corrections for 
E n k with k < n are zero. The form factors M n) „_2 receive additional corrections at order 0(q 2 ) from 
one-loop graphs with the insertion of pion operators, and E n ^ n receives corresponding contributions 
starting at order 0(q). 

For the form factors parameterizing moments of isoscalar GPDs, we find that B n ^, ^4 n ,fc and 
B Ut k at t = receive nonanalytic corrections of the form m 2 log(m 2 /// 2 ) from loops with nucleon 
operator insertions. No such corrections are found for A n ^ and C n . The form factors B n ^ n —2 

at t = 

receive in addition m 2 log(m 2 /// 2 ) corrections from loop graphs with pion operator insertions, and 
the corresponding nonanalytic contributions to C n give a term proportional to m. To leading chiral 
order, loop graphs with pion operator insertions are the only source of nonanalytic m 2 dependence 
for the derivatives of the form factors at t = 0. The derivative of M nj „_2 diverges like log(m 2 /// 2 ) in 
the chiral limit, and the derivative of C n like 1/m. 
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